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Abstract 

We report results from the MuLan measurement of the positive muon life- 
time. The experiment was conducted at the Paul Scherrer Institute using a 
time-structured surface muon beam and a segmented plastic scintillator ar- 
ray. Two different in-vacuum muon stopping targets were used: a ferromag- 
netic foil with a large internal magnetic field and a quartz crystal in a moder- 
ate external magnetic field. From a total of 1.6 x 10 12 decays, we obtained 
the muon lifetime = 2196 980.3(2.2) ps (1.0 ppm) and Fermi constant 
G F = 1.166 378 7(6) x 10" 5 GeV" 2 (0.5 ppm). 

The Fermi constant Gp is best determined by measurement of the positive muon 
lifetime r M . Following work by van Ritbergen and Stuart [1] and Pak and Czarnecki 
[2] — that reduced the theoretical uncertainty in the relation between the muon lifetime 
and the Fermi constant to 0.14 pm — two new experiments have been performed to 
measure jHH]. Herein we report the results from the part-per- million measurement 
of the positive muon lifetime by the MuLan Collaboration. 

1 Experimental setup 

The experiment was conducted at the Paul Scherrer Institute using a nearly 100% 
longitudinally polarized, 29 MeV/c, // + beam from the 7rE3 secondary beamline at 
the 590 MeV proton cyclotron. Incoming muons were stopped in solid targets and 
outgoing positrons were detected in a finely segmented, fast timing, scintillator array. 
The photomultiplier signals from scintillator detectors were recorded by 450 MHz 
sampling-rate, 8-bit resolution, waveform digitizers and read out by a high-speed 
data acquisition. 

One important feature of the experiment was the application of a time structure on 
the muon beam. The time structure consisted of a 5 /is-long, beam-on "accumulation 
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Figure 1: Cutaway drawing of the experimental setup: (1) vacuum beampipe, (2) 
stopping target, (3) halbach magnet, (4) scintillator array, and (5) beam monitor. 



period" followed by a 22 //s-long, beam-off "measurement period". The structure 
was imposed on the continuous beam using a custom-built, fast-switching, 25 kV 
electrostatic kicker. With the kicker high voltage off, the muons were transported 
to the target. With the kicker high voltage on, the muons were deflected into a 
collimator. The structure was important in avoiding the need to associate positrons 
with parent muons, a requirement that limited the statistics of earlier measurements. 

Another important feature of the experiment was the setup and the choice of the 
muon stopping targets. First, the 7rE3 beamline was extended through the positron 
detector with the stopping target mounted in the beam vacuum. This design reduced 
the number of upstream stops. Second, two combinations of stopping target materials 
and transverse magnetic fields were used: a magnetized Fe-Cr-Co foil (Arnokrome- 
III) with a 4 kG internal i?-field and a quartz crystal disk (SiO^) in a 130 G external 
B- field. In the ferromagnetic target, where muons reside as diamagnetic ions, the 
fi + precession frequency was ~50 MHz. In the quartz target, where muons reside 
mostly as paramagnetic atoms, the /i + e~ precession frequency was ~180 MHz. The 
fast precession yielded a roughly 1000-fold reduction in the ensemble-averaged fi + 
polarization via the spin dephasing during the beam accumulation. 

The positron detector was constructed of 170 triangle-shaped plastic scintillator 
pairs arranged in a soccer ball (truncated icosahedral) geometry (Fig. [T]). Each pair 
comprised an inner plastic scintillator tile and outer plastic scintillator tile. The pairs 
were grouped into ten pentagonal enclosures containing five tile-pairs and twenty 
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hexagonal enclosures containing six tile-pairs, which together formed the soccer ball 
geometry. The segmentation was important in reducing positron pile-up in individual 
detector elements. The symmetric arrangement of detector elements was important 
in reducing the effects of muon spin rotation / relaxation (/xSR). 

2 Data analysis 

A total of 1.0 xlO 12 decays from stops in Arnokrome-III were collected in our 2006 
production run and 0.7xl0 12 decays from stops in quartz were collected in our 2007 
production run. The runs yielded about 130 terabytes of digitized pulses. The anal- 
ysis was conducted using the 88-teraflop ABE cluster at the National Center for 
Supercomputing Applications (NCSA). 

The time and amplitude of pulses were determined from least square fits to the 
individual digitized waveforms. The procedure involved fitting a relatively high- 
resolution standard waveform (0.022 ns sampling-interval) to a relatively low-resolution 
individual waveform (2.2 ns sampling-interval). The high-resolution standard wave- 
forms were constructed by combining a large sample of 450 MHz sampling-rate, single- 
pulse, digitized waveforms. If the quality of the fit to a single pulse was not acceptable, 
the algorithm attempted to improve the fit by adding pulses or removing pulses. 

Two cuts were applied on the amplitudes and the times of the hits before con- 
structing inner-outer tile-pair coincidences. One cut defined an unambiguous software 
amplitude threshold A t h r for detector hits. Another cut defined an unambiguous arti- 
ficial deadtime ADT between detector hits. Hits that survived these cuts were sorted 
into time distributions of inner singles, outer singles and inner-outer coincidences. 
The construction of coincidence histograms with different thresholds A t h r and artifi- 
cial deadtimes ADT was important for studying such effects as pulse pileup and gain 
variations. 

2.1 Systematics effects 

If a hit occurs in the artificial deadtime of an earlier hit it is lost. Because these 
losses are more frequent during the high rates at early measurement times than the 
low rates at late measurement times, such pile-up distorts the time spectrum. 

Our procedure for correcting for pileup took advantage of the time structure of the 
incident beam. The pileup losses were statistically recovered by replacing the lost hits 
in each measurement period with measured hits at equivalent times in neighboring 
measurement periods. For example, to correct for leading-order pileup, if a hit is 
observed at time ti in fill j (the "trigger" hit), a hit is searched for within the interval 
U — >• U + ADT in fill j + 1 (the "shadow" hit). Adding the resulting histogram of 
shadow hit times to the original histogram of trigger hit times thereby statistically 
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recovers the lost hits. Similar methods were used to correct for higher-order pileup 
and accidental coincidences. 

In filling histograms only hits with amplitudes that exceed the threshold A t h r were 
used. Consequently, if the detector gain varies over the measurement period, then the 
time histogram will be distorted, either by additional hits rising above the amplitude 
cut or by additional hits falling below the amplitude cut. 

Our procedure for correcting for gain variations took advantage of the measure- 
ment of the positron pulse amplitude during the measurement period. The detector 
gain versus measurement time was determined from the MIP (minimally ionizing par- 
ticle) peak of the decay positrons in the scintillator tiles. Using the MlP-peak position 
versus measurement time, the time histograms were corrected for gain variations. 

2.2 Lifetime fits 

A simple procedure was used to extract the lifetime r M from the Arnokrome-III target. 
The summed tile-pair time histogram of coincidence hits was fit to Mr*/ T " + C. 
The approach relied on sufficient cancellation of Arnokrome-III //SR effects by spin 
dephasing and detector geometry. No evidence of /xSR distortion were observed in 
the time spectra from the Arnokrome-III target. 

A more detailed procedure was needed to extract the lifetime from the quartz 
target. First, geometry-dependent effective lifetimes were extracted for each tile-pair 
from fits to 

N(t) = Ne~ t ^ s [l + f(t)} + C, (1) 

where f(t) accounts for time-dependent effects of transverse- field (TF) spin precession 
/ relaxation. Then, the lifetime r M is extracted from the 170 effective lifetimes r e ff 
from a fit to 

T tS {.e B ,<f> B ) = T tt (i + s{e Bt <i> B )), (2) 

where 5(9 B ,<p B ) accounts for geometry-dependent effects of longitudinal-field (LF) 
spin relaxation. Together these steps accounted for all observed features of quartz 
TF/LF //SR. 

Note a number of datasets were accumulated with different magnetic field orien- 
tations and different positron detector offsets. These datasets changed the magnitude 
and orientation of TF/LF /iSR effects. The agreement between lifetime results for 
different configurations was an important verification of the fitting procedures. 

3 Results 

The individual results [5] for the muon lifetime from the Arnokrome-III dataset and 
the quartz dataset are 

r M (Arnokrome - III) = 2 196 979.9 ± 2.5(stat) ± 0.9{syst) ps (3) 
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and 

r M (quartz) = 2 196 981.2 ± 3.7(stat) ± 0.9(syst) ps. (4) 

The combined result 

r M (MuLan) = 2 196 980.3 ± 2.1(stat) ± 0.7{syst) ps (5) 

is obtained from the weighted average of the individual values with the appropriate 
accounting for the correlated uncertainties. The result corresponds to an overall 
uncertainty in the muon lifetime of 2.2 ps or 1.0 ppm — a thirty-fold improvement over 
earlier generations of lifetime measurements. The largest contributions to systematic 
uncertainties are associated with the aforementioned pulse pileup, gain changes and 
/iSR effects, as well as the knowledge of the time independence of the beam extinction 
during the measurement period. 

The final result for r M is in reasonable agreement with the earlier work of Duclos 
et al, Balandin et al, Giovanetti et al. and Bardin et al. However, it is in marginal 
disagreement (2.9 a) with the commissioning result from the FAST collaboration 
experiment jlj. 

We use the relation obtained by van Ritbergen and Stuart (vRS) pQ for the 
determination of the Fermi constant Gf from the measurement of the muon life- 
time t m . The vRS relation was derived using the V-A current-current Fermi inter- 
action with QED corrections evaluated to 2-loop order. It yields (jp(MuLan) = 
1.166 378 7(6) x 10~ 5 GeV~ 2 (0.5 ppm), a thirty-fold improvement over the 1998 Parti- 
cle Data Group [5] value that pre-dates the vRS theoretical work and recent lifetime 
measurements. The 0.5 ppm error is dominated by the 1.0 ppm uncertainty of the life- 
time measurement, with contributions of 0.08 ppm from the muon mass measurement 
and 0.14 ppm from the theoretical corrections. 
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